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ABSTRACT 


LOAD TEST ANALYSIS AT 
TOA INDUSTRIAL SITE, KAWASAKI, JAPAN 

by 

WILLIAM PATRICK DONNELLY 


Submitted to the Department of Civil Engineering on 19 August 
1968 in partial fulfillment of the requirements for the de¬ 
gree of Master of Science. 


Heavy loads on soft soil layers are becoming a common 
problem as good available sites are vised up. Much of the in¬ 
dustrial land around Tokyo has been reclaimed by pumping fill 
hydraulically over soft clay layers which underly Tokyo Bay. 
The weight of the fill and deep pumping of water in the area 
are causing area-wide compression of the clay layers. This 
settlement, combined with settlement under the weight of tanks 
on the reclaimed sites, is enough to cause serious problems. 

A test embankment was placed on the Toa Industrial Site 
400 in order to evaluate these settlements. Carrier (1966) 
attempted to predict before-hand, by one-dimensional analysis, 
the Skempton-3jerrum method and by the Lambe Stress Path 
method, the initial excess pore pressure and settlement due 
to this test load. The predictions were compared with actual 
field data and the differences rationalized by considering 
the methods of prediction. It was found that the stress path 
method gave the most reasonable results and that all methods 
of prediction improved when stress increments calculated by 
the finite element technique rather than elastic theory were 
used. 


Thesis Supervisor: T. William Lambe 

Title: Professor of Civil Engineering 
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CHAPTER I 


INTRODUCTION 

Heavy loads placed on soft soil layers constitute a 
problem area of soil mechanics. It is difficult to predict 
to what degree and how rapidly the load will settle. The 
purpose' of this thesis is to analyze an actual case illus¬ 
trating this problem by comparing predicted and measured 
performances. From this comparison, tentative conclusions 
are reached regarding the actual mechanism of settlement 
and true soil parameters. This comparison also allows ev¬ 
aluation of prediction techniques used for load settlement. 

Settlement under a load test on an industrial site in 
Kawasaki, Japan is the case to be analyzed. The soil pro- 
file of the site shows approximately 39 meters of soft clay 
overlaid by 11.4 meters of hydraulic fill. These two layers 
in turn rest on a rigid base of dense sand and stiff clay. 

An abundance of field data showing movements and pore pres¬ 
sures in each of the soil layers is available. 

Predictions of pore pressures and settlements under the 
load test were made by Carrier (1966) and by the author. 
Carrier used three methods .to predict pore pressures and 
settlements: One-Dimensional Analysis, Skempton-Bjerrum 

Method, and the Lambe Stress Path Method. The stresses and 
stress increments used in the soil tests associated with 
these methods were obtained from elastic theory. The author. 
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after studying the field data and comparing it with Carrier s 
predictions, used for his predictions the same three methods 
but with different soil parameters and stress increments ob¬ 
tained by the finite element technique. 

Comparisons were then made between the settlements and 
pore pressures, both predicted and actual, to determine which 
soil parameters and prediction technique best described the 
actual case. 
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CHAPTER II. 


DESCRIPTION OF THE TOA INDUSTRIAL SITE 
AND LOAD TEST 


A. Site History 

A severe shortage of industrial land in Japan has neces¬ 
sitated hydraulic reclamation of land, especially in the 
Tokyo area. The site under consideration in this thesis is 
owned by Toa Nenryo Kogyo K.K., situated in Kawasaki City, 
and was made by reclaiming land in Tokyo Bay. Figure II-l 
shows the location of the site (Site 400) which was filled 
during the period of April 1959 through January 1963 in two 
distinct phases. Using fill from both the canal and Bay, the 
northern and western sections were filled first, followed by 
the remaining portions. 

B. Soil Profile 

The soil profile of Site 400 is shown in Figure II-2. 

The base rock is mudstone overlaid with diluvial sediments 
extending upward to an elevation of about -50 meters. Desi- 
cation and overburden removal overconsolidated the diluvial 
soils. They are therefore relatively strong and incompres¬ 
sible, having little effect on the foundation engineering 
problems at the site. The only complication resulting from 
the diluvium concerns the sand-gravel layer which is an aqui¬ 
fer. Inland from Site 400, water for industrial use is pump¬ 
ed from deep wells extending down to this aquifer. 
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The soils extending from an elevation of approximately 
-50 meters up to the surface are alluvial sediments. The 39 
meters of soft clay have been divided into three layers on 
the basis of their engineering properties: Clay I, Clay II, 
and Clay III. The top stratum of Site 400 is the hydraulic 
fill, a clay-sand-silt mixture dredged from the Bay and canal. 

Study of the field data and tests run in the Soil Mech¬ 
anics laboratory at M.I.T. indicate that surface loads will 
affect only the fill and the three soft clays. Figure II-3 
shows a generalized soil profile of these soil strata and 
their properties as determined from the laboratory tests. It 
is emphasized that the laboratory test results for permeabi¬ 
lity and coefficient of consolidation of the fill are higher 
than those of the clay layers. It is expected that consoli¬ 
dation of the fill under a load would occur more quickly than 
would take place in the clay. 

C. Site Settlement 

The reclamation area including Site 400 is settling, 
without structural weight being imposed, for two reasons. 
Water is being pumped from the sand-gravel aquifer, causing 
excess pore pressure in the soft clay, which results in con¬ 
solidation of these layers as the pore pressure dissipates. 
The second reason is that consolidation of the same clay 
layers is still occurring under the weight of the placed 
fill. 
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Verification of these reasons can be obtained through 
examination of the piezometer data in the form of measured 
total head versus elevation as presented in Figure II-4. 

The pumping has caused the total head to drop below static 
head in the soil beneath Clay III while the addition of the 
fill has caused excess head above static in the clay layers. 

As these layers drain and consolidate, the total head will 
finally reach the line labelled steady, state pumping. Site 
settlement will be complete at that time. 

These two factors (consolidation due to pumping and 
consolidation under the fill) have resulted in 0.5 to 1.0 
meters of settlement to date around Site 400. Figure II-5a 
shows the whole area settlement and Figure II-5b gives the 
true rate of settlement at Site 400. 

Predictions have been made for the magnitude and rate of 
settlement caused by the two factors (Wong 1968, Lambe 1968) . 
It is estimated that more than three meters of settlement of 
the ground surface will take place due to consolidation of the 
clay layers; consolidation will not be complete until after 
1990. Comparisons of the predicted and actual rates of set¬ 
tlement are in very close agreement to date. 

D. The Load Test 

A test embankment was constructed on Site 400, as shown 
in Figure II-6, in order to gather data concerning bearing 
capacity and settlement as they relate to the subsoils in 
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that area. The theory was that the data could be used to 
predict the performance of structures (such as oil tanks) to 
be built on the site at a later date. 

The embankment was constructed in circular lifts of 
one meter in height to a final height of 6 meters and having 
a bottom diameter of 54.4 meters and top diameter of 40 
meters. Construction of the test embankment commenced on 
19 October 1966 and was completed 16 November 1966, each 
lift being placed over approximately four working days. Be¬ 
fore construction of the embankment began, special treatment 
was given to soil sections which were to underly the embank¬ 
ment. One third of the area was treated with 85 triangularly 
spaced sand drains, each 40 centimeters in diameter and ex¬ 
tending to a depth of 15 meters from the surface. In another 
third of the area, 235 paper drains were installed, also 
spaced triangularly and extending to a depth of 15 meters. 

No treatment was given to the remaining third portion of the 
area. 

In order to observe and monitor the effects of the test 
embankment on the underlying soil layers, instrumentation was 
installed. A total of 14 settlement plates, 10 piezometers, 
and 3 inclinometers were positioned as shown in Figure II-7. 
Settlement plates S-2, 3, 4, 5, 7, 8, 10, 11, 13, and 14 are 
screw type while the remaining settlement plates (S-l, 6, 9, 
and 12) are plate type. Piezometers P-6 and P-7 are the 
Casagrande type (open-standpipe porous tube). Piezometers 
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P-1, 2, 3, 4, 5, 8, 9, and 10 are open-standpipe piezometers 
of the Geonor type. Inclinometer 1-1 uses a stationary pen¬ 
dulum and electrical resistance. Inclinometers 1-2 and 1-3 
are the Wilson type, employing electrical resistance and a 
pendulum which slides up and down the pipe. 
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CHAPTER III 


PERFORMANCE OF THE LOAD TEST 

A. Horizontal Movement 

A great deal of instrumentation was installed under the 
test embankment in order to measure both vertical and hori¬ 
zontal movements of and pore pressures in the soil layers. 

The inclinometer readings measuring horizontal soil movement 
did not fluctuate from their initial values, indicating in¬ 
valid readings or no horizontal displacement, probably in¬ 
valid readings. 

B. Vertical Movement 

From the movement of the settlement plates it is pos¬ 
sible to determine the compression of each of the soil layers 
and to compare these compressions with those of the same 
layer in an area with different treatment. Figures Ill-la, 
b, c, and d show the compression of the soil layers in each 
of the different areas. It is seen that the compression of 
Clay I, II, and III in Figures Ill-la and Ill-Id differ mar¬ 
kedly from the compression of the same layers in Figures III- 

lb and III-lc. This difference is so great that it cannot 

♦ 

reasonably be attributed solely to actual compression differ¬ 
ences in the different sections. Based on soil properties and 
predictions from soil mechanics theory, it is concluded that 
the settlement curves showing compressions in the untreated 


13 







and paper drain sections shown in Figures Ill-lb and III-lc 
represent actual movements and that settlement plates S-3 and 
S-14 (under the center and in the sand drain section) are 
reading incorrectly. The reason for the suspected error in 
settlement of these plates is that they were probably carried 
down by friction. Therefore, settlement curves based 'on com¬ 
pression of soil layers in the center and sand drain sections 
measured by these plates are not used, except to illustrate 
the large differences. 

Figures III-2a, b, and c compare the compressions of 
the same soil layers in the paper drain and untreated sec¬ 
tions; also, it is seen how much difference is caused by the 
faulty settlement plates (S-3 and S-14) in the center and 
sand drained sections. These figures indicate that special 
draining treatment has little effect on the magnitude or 
rate of settlement in any of the soil.layers. 

C. Pore Pressures 

Piezometers were installed to measure pore pressures 
under the center and in each of the treated sections under the 
load test. (Piezometers P-5 and P-6 were either suspected 

damage during construction or ceased functioning on their own 

* 

and were replaced by piezometers P-5a and P-6a. Piezometer P- 
10 was also discontinued but was not replaced.) Figures III- 
3a and III-3b show the pore pressures (as measured by these 
piezometers) versus time from date test embankment was placed. 
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It can be seen from these plots that the piezometers in the 
clay have a time lag of two months or more (from the time of 
embankment placement) before pore pressure buildup is regi¬ 
stered. Also, piezometer P-5a has given very erratic read¬ 
ings which indicate it may be unreliable. 

D. Immediate Settlements 

Soft soil layers, when subjected to a heavy load, will 
undergo immediate settlement which is caused by a change in 
shape but not in volume of a soil element; there is no drain¬ 
age or consolidation involved. It is possible, in this case, 
to estimate the magnitude of immediate settlement or compres¬ 
sion of each of the soil layers by determining, from settle¬ 
ment plots such as Figures Ill-lb and III-lc, the compression 
of these layers at the middle of the loading period of the 
test embankment. Figures III-4a and III-4b show the same 
figures replotted to a larger time scale where 0 time is the 
middle of the loading period. 

The drawback to this method of estimating immediate 
settlements is that the choice of the middle of the loading 

period as the point of average immediate settlement is 

* 

somewhat arbitrary. In this case it is based on the assump- 

* 

tion that a load placed uniformly over a specific period of 
time can be closely approximated by the full load placed in¬ 
stantly at the middle of the period (Taylor 1948). A slightly 
different choice of the time of average immediate settlement 
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causes a large difference in magnitude of the estimates due 
to rapid compression of the soil layers during this time. 
However, in Figures III-4a and III-4b the curves have been 
extrapolated back to zero time rather than forward from the 
beginning of placement of the embankment. Therefore, the 
trend of the curves is due to consolidation and the immediate 
settlement, as shown by the intersection of these curves with 
zero time, should be more accurate. The results are shown in 
Table III-l. 

E. Initial Excess Pore Pressure 

From the piezometer data, it is possible in this case to 
make an estimate of initial excess pore pressures. Before 
this can be done, the pore pressure curves as shown in Fig¬ 
ures III-3a and III-3b (since there is a time lag) must be 
extrapolated back to zero time. Again, the zero time is 
chosen as the middle of the loading period. As a first esti¬ 
mate, straight lines are used in the extrapolation. The ini¬ 
tial excess pore pressures are the differences between the 
piezometer readings at the end of September 1966 and the ex¬ 
trapolated values at the middle of the loading period. They 
have been tabulated in Table III-2. 

F. Percent of Consolidation 

The pore pressure distribution due to the load test in 
the soil layers over a period of a year, as determined from 
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the straight line extrapolation in Figure III-3a, has been 
plotted in Figure III-5. Using these figures, it is possible 
to determine (as shown in Table III-3) the degree of consoli¬ 
dation for every three months throughout the first year. 

In this case, the extrapolation also has been up to the 
judgement of the individual. A slightly different angle of 
the straight line approximation can result in large differen¬ 
ces in initial excess pore pressure. 

One method of checking the reliability of the piezo¬ 
meter data and the accuracy of the straight line extrapola¬ 
tions to zero time is to compare the percentage of consolida¬ 
tion values during the year period to the actual settlement 
curves. For example, Clay I in the untreated section has com¬ 
pressed 57 centimeters in a year. If 57 centimeters less im¬ 
mediate settlement represents 69% of the consolidation under 
the test embankment as suggested by Figure III-5 and Table 
III-3, then the settlement curve based on the percent consoli¬ 
dation values for Clay I would appear as shown in Figure III- 
6a. Also shown in Figure III-6a is an actual curve for Clay I 
The same procedure has been carried out for Clay II and III 
and the results are shown in Figure III-6b. It is seen that 
consolidation as determined from the straight line extrapola¬ 
tion does not occur as rapidly as indicated by the actual 
curve. The most probable reason is that a straight line ex¬ 
trapolation is not a good approximation of the actual case. 

The extrapolations are probably more like those shown in 
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Figures III-7a and III-7b. (It was decided from Figure III-5 

that the data from piezometer P~5a was too high and that P-5 
gave a better representation of the actual case.) Using these 
figures to find new values for percent consolidation and ap¬ 
plying them to a settlement curve results in the third curve 
shown in Figure III-6a and III-6b. Clearly, this curve is 
much closer.to the actual settlement curve. Therefore, the 
curved line extrapolations in Figures III~7a and III-7b are 
closer to reality. Tables III-4 and III-5 give new values of 
initial excess pore pressures and percent consolidation based 
on the revised extrapolations in Figures III-7a and III-7b. 
Figure III-8 shows the revised pore pressure distribution 
under the center of the load test. 


18 




CHAPTER IV- 


COMPARISON OF PREDICTED AND'MEASURED 
PERFORMANCE OF THE LOAD TEST 

A. Introduction 

Carrier (1966) predicted settlements and pore pressures 
along the center line under the test embankment by four 
methods: field data from embankments on Site 200, one¬ 

dimensional analysis, the Skempton-Bjerrum method, and the 
Lambe Stress Path method. He assumed that excess pore pres¬ 
sures comparable to those in the clays were also built up 
in the fill layer (undrained loading) and used this fact in 
his predictions. However, it was seen in Chapter III and 
Figure III-8 that there was no registered pore pressure build¬ 
up in the fill and therefore a drained loading probably more 
closely approximated the actual case. For these reasons the 
author directed most of his attention to a comparison of pre¬ 
dicted and measured performance in the clay layers. Carrier's 
predictions, by the three theoretical methods, of amount and 
rate of compression of the soil layers are compared with the 
actual compressions in Figures IV-la, b, and c. Shown in 
Figure IV-2 is the comparison of Carrier's predictions of ini¬ 
tial excess pore pressure and the author's estimation of the 
actual initial excess pore pressures (Chapter III). Finally, 
Table IV-1 is a comparison of the predicted and measured im¬ 
mediate settlements. 
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B. Discussion of the Results of Comparison 

1. Rate of Compression 

The most obvious difference between the predicted 
and measured values of compression is the rate of compression 
in the clay layers. The predicted rates were based on a 
laboratory coefficient of consolidation c v determined under 
one-dimensional conditions. Since the predicted rates were 
much slower than the measured rates of compression, it is 
certain that radial drainage plays a large role and the one¬ 
dimensional role does not come close to approximating the 
actual consolidation. 

Unfortunately, the theory of three-dimensional consoli¬ 
dation is not satisfactorily advanced to be of use in this 
case. One method of overcoming this problem is to use a 
higher value of c v than obtained in the one-dimensional case. 
Wong (1968) backfigured from the load test a c v approxi¬ 
mately 8 times greater than the c v obtained in the labora¬ 
tory. Further proof that the settlement under the test em¬ 
bankment is not one-dimensional is the fact that the rate of 
area-wide settlement, which is likely to be one-dimensional, 
compares closely with the rate predicted on the basis of one¬ 
dimensional consolidation. 

2. Magnitude of Compression 

Since the test embankment has been in place less 
than two years, and primary consolidation is probably not yet 
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complete, an extrapolation must be used to estimate the final 
amount of compression of each clay layer. This was done by 
determining the percentage dissapation of excess pore pressure 
at a certain time from Figures III-7a, III-7b and III-8 as 
listed in Table III-5. With the percentage dissapation and 
measured amount of compression at the same time, it is- pos¬ 
sible to estimate the amount of compression at 100% dissapa¬ 
tion of excess pore pressure. This has been done for each 
clay layer and the estimated values are compared with Car¬ 
rier's predicted values in Table IV-2. 

In order to establish reasons for the differences be- 
tween the predicted and estimated final values of compres¬ 
sion, the methods of prediction must be examined in more de¬ 
tail. 


In Carrier's one-dimensional analysis, he used the fol¬ 
lowing formula for primary consolidation: 


p c “ ! lfe - lo 9 dH 


o 


VO 


To determine values of ° v f ' a< ^ded a vertical stress in¬ 
crement Aa = 10.93 t/m 2 to o „ for each clay layer. This 
v VO 

stress increment is greater than those determined by elastic 
theory and by the finite element technique, especially in the 
lower layers of Clay II and III. A graphical representation 
of the vertical stress increments determined by these three 
methods is given in Figure IV-3. Certainly the stress incre¬ 
ments determined by elastic theory and the finite element 
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technique are more realistic; this .is supported by the com¬ 
parison of predicted and estimated values of compression. 

The compression predicted for Clay I by one-dimensional means 
is slightly higher than the estimated value, while the com¬ 
pression predicted for Clay II and III is almost twice as 
great. The much larger value predicted for Clay II and III 
is probably largely due to the greater stress increments pre¬ 
dicted by this method and used in consolidation tests on 
samples from these layers. 

. The Skempton-Bjerrum method (1957), developed for eval¬ 
uating settlement in clays, is a semi-empirical method but 
is more realistic than one-dimensional analysis. This method 
recognizes that in settlement problems in clay, the total set¬ 
tlement has three components: 

p i + p c + p s 
initial settlement 

primary consolidation settlement 

secondary consolidation settlement 

The initial (immediate) settlement is caused by a soil ele¬ 
ment changing shape but remaining at constant volume; the soil 
element undergoes undrained shear. Usually, is computed 

by the theory of elasticity, assuming the soil mass to be iso¬ 
tropic and homogeneous. However, Carrier used a different 
method. He divided the soil profile into 26 layers, collected 
all the triaxial data available for these layers, normalized 


where 


P 

P. 
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the data for each test, and prepared an average plot of strain 
from the undrained portions of the tests for each clay-. He 
then determined the immediate settlement by multiplying each 
strain by the thickness of its corresponding layer. The re¬ 
sults were shown in Table IV-1 and are good if the amount of 
settlement is considered. 

For primary consolidation settlement, Skempton and Bjer- 
rum used the same expression as for one-dimensional consoli¬ 
dation with the following modifications: 


= f 


or 


where 


m 


v 


where 


so that 


1+e 


log 


’vf 


0 vo 


dH 


f m v Aa v dH 


lQ g a vf /a vi 


1+e 

o 


A a 


v 


and substituting Au for A o 


v 


Au = initial excess pressure 
= Au h + A(Aa v “Aa H ) 


/ m Au dH 
v 


This recognizes that settlements are controlled by Aa H as 
well as Aa v . 

The pore pressure parameter A is usually selected from 
a chart which gives a constant range of values of A for dif¬ 
ferent types of clays. However, it has been shown that A is 
not constant for a given soil but depends on strain, amount of 
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over-consolidation, anisotropic consolidation, type of shear, 
and environment. 

Therefore, Carrier used a different method to select 
values of A . He calculated values for each of the original 
26 layers from the normalized triaxial data. A comparison of 
his predictions of Au for the Skempton-Bjerrum method and 


the measured values was shown in Figure IV-2 and will be dis¬ 
cussed in the next section. It is seen that he under-predicted 
values of Au for Clay I. If Clay II and Clay III are taken 
as one layer, and the predicted and measured values averaged, 
his prediction for Clay II and III is reasonably close to the 
measured value. 

After determining u on a three-dimensional basis, the 
Skempton-Bjerrum method assumes the consolidation settlements 
are one-dimensional and uses the results of oedometer tests 
to determine the amount. 

Thus: 


m 


v 


1+e 


log a c /o . 
u vf vi 

Aa 


where 


A a = a £ - a . 

v vf vi 


and a . is the vertical stress acting on the soil after un- 
vi 

drained shear has occurred (Figure IV-4). It is not the same 

as a , the in situ vertical stress used in the one-dimens- 
vo 


ional analysis 
Since a 


, - o • = Au 

vf VI 


then 


A a v = Au 
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and the two cancel out in the formula for 


thus: 


= / 


1+e 


log dH 


Carrier used his predicted value of Au for Clay I, 
which was lower than the field value as estimated by the 
author, in determining for the preceding formula. 

This helps to explain why, using the Skempton-Bjerrum me¬ 
thod, he predicted such a low value of compression for Clay 
I as compared to the field value estimated by the author 
(Figure IV-2). 

Both the previous methods discussed, the one-dimensional 
analysis and the Skempton-Bjerrum method, have the same draw¬ 
back: neither method really approximates what happens in the 

field. Figure IV-4 shows a first approximation of the actual 
stress path that is followed by a soil element under the load 
test. It also shows the stress paths that were followed in 
laboratory tests on soil samples for both methods of predic¬ 
tion discussed so far. 

If the strains between the same vertical stresses were 
the same for different stress paths, the Skempton-Bjerrum 
method would be a simple and accurate way of determining 
settlements in clays if au could be predicted accurately. 
Unfortunately, the strains for two different stress paths be¬ 
tween the same vertical stresses are not usually the same, 
as shown by Carrier (1966) in Figure IV-5. 

The fact that different stress paths can result in large 
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differences in strain, as seen in Figure IV-5, is another 
factor which contributed to the differences between actual 
and predicted values of magnitudes of.compression. 

In some cases, however, it is possible that this factor 

can compensate for the other errors. For example, in the 

Skempton-Bjerrum method where the Au and thus a are 

under-predicted for Clay I, leading to a prediction of a 

smaller p , the larger strain given by the sample in the 
c 

oedometer during the consolidation portion will lead to a 
larger prediction of strain than if the sample were run along 
the actual stress path. 

To overcome these difficulties, Lambe (1967) developed 
the Stress Path method of predicting settlements where p c 
is not based on oedometer tests but on. three-dimensional tri- 
axial tests which more accurately duplicate field conditions. 
An "average" soil sample is first consolidated to in situ 
stress conditions then loaded under the same stresses ex¬ 
pected in the field at that depth (the stress increments are 
usually determined from elastic theory)--first in undrained 
shear, then consolidated at constant shear stress (Figure IV- 
4). The strain is then measured and multiplied by the thick¬ 
ness of the soil layer to determine settlement or compression. 

In this case. Carrier determined initial excess pore 
pressure and immediate settlements in the same manner for both 
the Skempton-Bjerrum and Lambe methods. 

Carrier selected the average soil samples, stresses, and 
stress increments on the basis of a 60-meter diameter 
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embankment. After he completed his tests, the diameter of . 
the proposed embankment was changed to 40 meters, making it 
necessary to normalize his data and assume the vertical strain 
was proportional to log p^/p^ • 

It is seen from Table IV-2 that the Lambe method gives 
prediction of magnitude of compression of the soil layers 
closest to the field values. Carrier (1966) lists the test¬ 
ing difficulties for this method which undoubtedly had some 
effect on the results. 

3. Initial Excess Pore Pressure 

Figure IV-2 compared Carrier's predicted values of 
initial excess pore pressure with those estimated by the 
author from field data. As was stated earlier and seen in 
Figure IV-2, the predicted value of excess pore pressure for 
Clay I is considerably lower than the field estimated value. 
There could be several reasons for this. The possibility of 
faulty piezometer data should not be ruled out, although the 
data for P-4 in Clay I is relatively consistent and the ex- 

r 

trapolation to determine the initial excess pore pressure is 
supported by the field data for compression of Clay I. More 
likely, the error comes from the samples which were used in 
the triaxial tests upon which Carrier based his pore pressure 
predictions. Clay I appears to contain many small isolated 
pockets of sandy silt. When Clay I is taken as a whole, the 
pockets have little effect in decreasing pore pressure build¬ 
up or increasing dissapation. However, these pockets in a 
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small triaxial sample do have the effect of decreasing initial 
excess pore pressure. Therefore, it is probable that initial 
pore pressures measured in the triaxial tests were low. 

Another factor which can cause differences may be the 
fact that the stress increments in Clay I, due to the load 
test, bring Clay I under the center line very close, if not 
all the way, to failure (Chapter VI). If this happens in the 
field, the shear stresses are re-distributed along the center 
line, stopping further strain, and the pore pressures can con¬ 
tinue to rise. In a triaxial sample reaching failure under 
undrained loading, pore pressures will stop rising and fall 
as the strain increases. Therefore, pore pressures obtained 
in the field in this way cannot be measured in the laboratory. 
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CHAPTER V 


STRESS CALCULATIONS 


A. Introduction 

Ali of the methods which were used to predict settle¬ 
ments require that the stress increments due to the load test 
be estimated at various depths below ground surface. In the 
past, the most common way to do this has been by means of 
elastic theory in which the soil is assumed to be a homo¬ 
geneous, isotropic, elastic halfspace. This assumption is 
invalid for most soils, thereby introducing errors into the 
calculation of stress increments, especially in the horizon¬ 
tal direction. 

In the author's opinion, a better approach (which has 
been developed for use with computers) is the finite element 
technique because it can handle stratified soils whose layers 
have different moduli and Poisson's ratios. It therefore ap¬ 
pears to be a much closer approximation to the actual case. 

B. Moduli Used 

Both elastic theory and the finite element technique re¬ 
quire that Young's moduli and Poisson's ratios be determined 
in order to calculate the stress increments. Since elastic 
theory assumes that the soil is a homogeneous, isotropic, 
elastic halfspace, only a single value of Young's modulus and 
Poisson's ratio can be used. 
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For the finite element technique, a Young's modulus and 
Poisson's ratio is required for each different material used 
in the finite element analysis. The most significant errors 
which can occur using this technique are usually the result 
of incorrect moduli. 

According to Ladd (1964, 1965), undrained moduli are de¬ 
pendent on several factors including stress level, anisotropy, 
and overconsolidation ratio. Wong (1968) discusses methods 
of obtaining undrained moduli for the soil layers under the 
load test considering these factors. He calculated undrained 
moduli by determining stress level and then using the finite 
element technique. He also backfigured them, using the im¬ 
mediate settlements of the load test and triaxial data on the 
Kawasaki clays reported by Lambe (1962). 

In the field, the layer of fill drained much faster under 
the load test than the clay layers, as shown by the pore pres¬ 
sure data in the preceding chapters. Therefore, all shear 
which took place in the fill was considered drained shear and 
an undrained modulus was not determined for this layer. 

Based on this assumption, Wong assigned to the fill a 
Poisson's ratio of 0.33 and backfigured from the load test a 
Young's modulus of 150 t/m^. 

Wong also obtained drained moduli for each of the clay 
layers using elastic theory and assuming one-dimensional 
loading. 

Tables V-la and V-lb summarize the values of moduli and 
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Poisson's ratio obtained by using the various techniques on 
the Kawasaki soils. 

C. Stress Increments 

In order to use elastic theory to determine stress in¬ 
crements, the soil layers must be lumped together so that a 
single average modulus and Poisson's ratio may be selected. 
Also, the test embankment must be replaced by an equivalent 
cylinder (Figure V-l). This was done for different values of 
Poisson's ratio. The stress distributions which result are 
seen in Figure V-l, showing that the vertical stress incre¬ 
ment is independent of v while the horizontal stress in¬ 
crement is not. 

As stated before, the advantage of the finite element 
technique is its ability to handle problems involving soil 
layers having different moduli and Poisson's ratios. The 
input data for the computer is in the form of moduli and 
Poisson's ratio for each different material and numbers de¬ 
scribing elements and nodal points of a mesh superimposed -on 
the cross-section of the soil layers. The mesh used for this 
particular problem is shown in Figure V-2. A computer pro¬ 
gram called FEAST I was used which was developed by E. L. 
Wilson of the University of California in 1966 and modified 
by John Christian of M.I.T. The program allows close simula¬ 
tion of boundary conditions, such as the rigid layers under 
the soft clays. 

The distribution of stress increments obtained by this 
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method under undrained and drained conditions are shown in 
Figures V-3a and V-3b, respectively. Figure V-4 compares the 
elastic stress distribution with stress distributions obtained 
from the finite element technique, using both undrained and 
drained moduli and corresponding Poisson's ratio. Some im¬ 
portant points to notice from this figure are the small dif¬ 
ferences in AcJi , as determined by the previously described 
methods, and the relatively large variance in Ac 3 . An er¬ 
ror in calculating Act 3 will have no effect on the predic¬ 
tion of settlement by one-dimensional analysis. However, it 
will affect the determination of Au needed in the Skempton- 
Bjerrum method of calculating settlements. The Lambe Stress 
Path method will be the most affected, however, since the 
terminating point of the stress path along which laboratory 
samples are taken in a triaxial test depends on both Aa 3 
and Aa 1 . 
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CHAPTER VI 


REVISED PREDICTIONS OF INITIAL EXCESS PORE PRESSURES 
AND COMPRESSIONS OF SOIL LAYERS 

A. Introduction 

Analysis of the field data from the load test and com¬ 
parison of these data with previously predicted performance 
illustrate what probably happens in the subsoil. Consid¬ 
ering this, it is possible to revise certain aspects of the 
methods of prediction with new techniques (such as finite 
elements) in order to make better predictions. Suggestions 
are made for improvements of each of the three methods dis¬ 
cussed previously. However, most of the effort has been spent 
on the Stress Path method since it is considered more real¬ 
istic and straight-forward than the others and should even¬ 
tually yield the best predictions. v 

B. ' One-Dimensional Analysis 

If the diameter of the load test were large (relative 
to the thickness of the compressible strata), the one¬ 
dimensional method of prediction would be a good and simple 
approximation of the field case. This is not true in this 
example because consolidation has been shown to be three- 
dimensional. An improvement over Carrier's predictions, how¬ 
ever, could be made in determining Aa v by elastic theory or 
the finite element technique, rather than allowing Aa v , at 
all depths, to be equal to the stress increment added at the 
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surface. This has been done for the load test with the re¬ 
sults summarized in Table VI-1. These predictions are seen 
to be closer to the actual case than Carrier's original esti¬ 
mate. (Compare pore pressures with those in Table III-4 and 
settlements with those in Table IV-2.) 

C. Skempton-Bjerrum Method 

The Skempton-Bjerrum method of predicting settlements, 
as stated in Chapter III, is a better approximation of actual 
field conditions in this case. It considers immediate and 
secondary as well as primary compression, and comes closer to 
accounting for three-dimensional effects by using Aa 3 in 
determining Au and by using Au in calculating settlements 

Au = Act 3 + A (Act j -Act 3 ) 

p = H m Au 

K c . v 

Skempton and Bjerrum recommend that parameter A be se¬ 
lected from a chart which, in the author's opinion, is very 
limited (especially for the inexperienced) since it neglects 
specifically the'factors mentioned in Chapter III to which A 
is sensitive, and allows them only to be included intuitively 
The author agrees with Carrier (1966) and Wong (1968) that 
the parameter A be determined from previous CIU and CAU 
tests on clay samples taken from Site 200. Values of A de¬ 
termined by this method are shown in Table VI-2. 

Figure V-4 illustrated the relatively large differences 
in the Clay I layer between the Aa 3 determined by elastic 
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theory and the Aa 3 determined by the finite element tech¬ 
nique using undrained properties. The finite element tech¬ 
nique gives higher values of Aa 3 in Clay I, resulting in 
higher predicted values of An for that layer and bringing 
the predictions closer to the field measurement of Au . 

This supports the author's belief that the finite element 
technique yields a better representation of stress distri¬ 
bution than does elastic theory. Use of a higher value of 
Au in the one-dimensional formula employed by Skempton and 
Bjerrum for consolidation settlements should also give a 
larger value of compression in Clay I, bringing this predic¬ 
tion closer to the field data. A summary of the revised 
Skempton-Bjerrum predictions is given in Table VI-3. 

The prediction of Au and compression in Clay I is 
still lower than the field values. This may indicate that 
the A factor calculated for Clay I from triaxial tests on 
samples from Site 200 is not close to the A factor for 
Clay I in Site 400. 

However, the predicted values of initial excess pore' 
pressure, au , and compression for Clays II and III are re¬ 
latively close especially using the A factor determined 
from CAU tests. 

D. Lambe Stress Path Method 

The stress path concept of predicting settlements and 
initial excess pore pressures uses tests which better 
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duplicate field conditions than those used in the other meth¬ 
ods.’ On this basis it is expected that the stress path meth¬ 
od would yield better predictions. As was seen in Chapter IV, 
this was the case in Carrier's predictions. However, for the 
stress path method, stress increments determined by elastic 
theory were again used and, as noted previously, there are 
relatively large differences between horizontal stress incre¬ 
ments calculated by the finite element technique. The stress 
path method will be the most affected by changes in A a 3 
since the stress path along which samples are tested depends 
on Aa 3 as well as Ag x . 

Based on the assumption used in improving the previous 
prediction techniques (that the stress distribution deter¬ 
mined by the finite element technique is more accurate than 
that provided by elastic theory), new effective stress paths 
for each of the clay layers were drawn. These are shown in 
Figures Vl-la, b, and c. The end points, E , were deter- 
mined from the undrained stress increments AGi and Aa 3 , 
as shown in Figure V-4 and listed in Table VI-4. The undrained 
portions of the paths in Figure VI-1, E q -E^ ' were determined 
by estimating immediate strain for each layer under the load 
test, entering Carrier's curves of q/a„_ versus axial strain 
at the appropriate strain to determine q , then from his nor¬ 
malized curves of p versus q , determining p for that q 
(Table VI-5). Carrier showed several normalized curves of q 
versus axial strain and p versus q for each clay layer. 

For the preceding manipulations, it was therefore necessary 
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to determine from his data average curves for each layer; 
these are shown in Figures VI-2a, b, and c. 

Carrier assumed that vertical consolidation strain at 

constant shear stress in a triaxial test could be normalized 

similarly to vertical strain in an oedometer test. Based on 

this assumption, he plotted curves of p^/p^ versus vertical 

strain from Kawasaki Site 400 clay samples. Using this data 

it is then possible to calculate strains for consolidation at 

constant shear stress as takes place in the proposed stress 

paths. The predictions of p and Au based on these stress 

c 

paths are given in Table VI-6. It is noticed that these pre- 
dictions, when compared with the field data in Tables III-4 
and IV-2, are closer than the predictions made by Carrier 
based on stress paths using stress increments determined 
from elastic theory. 

An example of a field total stress path for the effec¬ 
tive stress paths just constructed would be as shown in 
Figure VI-3 for Clay I, where both Aa x and Acr 3 increase 
during undrained shear and then remain constant during con¬ 
solidation . 

During the undrained portion of the stress path, the 
soil element (sample) changes shape but not volume (see Fig¬ 
ure VI-4a). It is quite reasonable, therefore, to assume a 
Poisson's ratio of 0.5. However, during consolidation (the 
drained part of the stress path), the soil element changes 
shape and volume (see Figure VI-4b). Poisson's ratio decreases 
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during consolidation to a value of less than 0.5; this causes 
a change in total stresses. These facts are reflected in the 
stress increments determined by the finite element technique 
(see Figure V-4). Assuming a Poisson's ratio of 0.33 gives 
a smaller value of Aa 3 than for a Poisson's ratio of 0.5, 
especially in Clay I. Poisson's ratio, however, has compara¬ 
tively little effect on the vertical stresses. The change, 
then, in vertical effective stress during consolidation is 
equal to the initial excess pore pressure while horizontal 
effective stress during consolidation decreases with the 
total horizontal stress. Probably a better representation 
of the effective and total stress paths for an element in 
the middle of the Clay I layer is as shown in Figure VI-5. 

The stress values are given in Table VI-7. 

To test this theory, the author attempted to duplicate 

this stress path in triaxial tests on.samples from near the 

middle of Clay I. The undisturbed samples came from boring 

1-16, located as shown in Figure II-6. The stress path tests 

were conducted according to the steps outlined by Guertin 

(1967) except that the stress-controlled K q consolidation 

was carried out in a.strain control loading frame so that 

the triaxial cell did not have to be moved between the stress- 

controlled and strain-controlled portions of the test. This 

modification eliminates disturbance to the sample caused by 

the move. The samples were consolidated anisotropically along 

the K line to three times the in situ vertical effective 
o 

stress in order to eliminate the disturbance effects of 
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sampling procedures and handling. 

It was expected that effective stress paths for the un¬ 
drained portion of the tests would roughly approximate those 
of Carrier's, as shown in Figures yi-3 and VI-5 ( E^-E^ ). 

This, however, was not the case. Figure VI-6 ( E q -E^ ) shows 
the undrained portion of the stress paths as recorded in 
these tests. There are two reasons which may account for the 
difference. First, as mentioned earlier, the layer of Clay I 
contains many small isolated pockets of sand. It is possible 
for a triaxial sample to consist of as much as 50% sand which 
is certainly not representative of Clay I. This may have 
been the case with Carrier's samples, causing his results to 
differ from the author's. The other reason for the difference 
may be the improvement of triaxial equipment in the M.I.T. 
laboratory between the time of Carrier's tests (1965-66) and 
the author's. The author feels that his results better ap¬ 
proximate what actually happened in the field. 

The results of these tests were normalized and are shown 
in Figures VI-6, VI-7 and VI-8, and summarized in Table VI-8. 

The initial excess pore pressure obtained in these tests 
were much higher than expected from Carrier's data but very 
close to the estimated field value for Clay I. Strains in 
the samples tended to go slightly higher than would be expec¬ 
ted on the basis of field measurement. This could be explained 
by the fact that the stress path comes very close to the fail¬ 
ure envelope causing greater strains in a triaxial sample than 
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in the field where the stresses probably redistributed, there 
by retarding further strain. 

It is felt that this stress path is more nearly correct 
than those previously discussed because it does consider the 
change in Poisson's ratio between drained and undrained shear 
However, this type of stress path test must be carefully con¬ 
trolled, especially near failure, if it is to give accurate 
results. 
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CHAPTER VII 


CONCLUSIONS 


A. Introduction 

Table VII-1 gives a summary of the predictions of initial 
excess pore pressures and compression due to primary consoli¬ 
dation in each of the soil layers. The same quantities de¬ 
termined from the field data are also shown. 

It is seen that for each method, use of the finite element 
technique in determining stress increments produced an im¬ 
provement in the predictions as compared with the field values. 
The author believes that the finite element technique is su¬ 
perior to elastic theory in calculating stress increments. 

This table also shows that in using both elastic theory 
and the finite element technique, the Lambe Stress Path method 
gives predictions closest to the field values. The author 
feels that the stress path method is best because it is the 
only one of the three which attempts to duplicate field condi¬ 
tions. For this reason any improvements in techniques, such 
as in calculating stress increments, are adapted with more 
significance into the stress path method than into the others. 
As improvements such as these are added, predictions by the 
stress path method will converge on the true answer, while 
incorporation into the other methods may or may not. Certain 
aspects in the field, such as greater strains when a soil ele¬ 
ment nears the failure envelope, will not be reflected in 
tests used in the one-dimensional and Skempton-Bjerrum 


41 ’ - 






methods but will be included in the Stress Path method. 

B. Areas of Greatest Uncertainty 

As is evident throughout this thesis, the greatest uncer¬ 
tainty in predicting pore pressures and settlements under 
heavy loads in soft clay layers lies in the determination of 
stress increments in the subsoils. Although the finite ele¬ 
ment technique appears to be an improvement over elastic the¬ 
ory in this area, moduli (both drained and undrained) and 
Poisson's ratios are needed for each soil layer. Methods of 
determining these quantities are subject to error. 

A danger in all soil tests, especially the stress path 
tests, is not choosing a representative soil sample. There¬ 
fore, care must be taken to insure that the soil profile is 
accurate (concerning relatively homogeneous layers) and that 
the samples to be tested are truly representative of these 
layers. 

C. Suggestions for Future Work 

1. Field Load Tests 

When running a load test for the purpose of evalu¬ 
ating settlements and pore pressures, it is evident that great 
care must be taken in installing instruments (piezometers, 
settlement plates, and inclinometers) to give reliable data. 
Data from the inclinometers would have been valuable in eval¬ 
uating horizontal stresses under the load. Even more 
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valuable would be a load cell'which, could measure both ver¬ 
tical and horizontal total stress in the field without in¬ 
fluencing them too much. 

2. Laboratory Tests 

The author believes that using a combination of 
drained and undrained stress increments to determine the 
stress paths, as was described in the final part of Chapter 
VI, gives a more accurate representation of field conditions 
and thus yields better predictions. The stress path tests 
conducted in this manner gave good results but more tests 
are needed to confirm this. 

It is seen from Table VII-1 that all methods of predic¬ 
tions give reasonable results. For this reason it is felt 
that for simple, inexpensive jobs where the amounts of set¬ 
tlement are not critical, one-dimensional analysis and the 
Skempton-Bjerrum method may be preferred to the Stress Path 
method since stress path tests are more time consuming and 
expensive to run. 
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